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Abstract. We study the tidal evolution of a binary system con- 
sisting of a 1.4 Mq compact object in elliptic orbit about a 
10 Mq uniformly rotating main sequence (MS) star for vari- 
ous values of the initial orbital parameters. We apply our pre- 
viously published results of 2D non-adiabatic calculations of 
the non-radial g- and r-mode oscillations of the uniformly ro- 
tating MS star, and include the effects of resonant excitation 
of these modes in the tidal evolution calculations. A high or- 
bital eccentricity enhances the effectiveness of the tidal inter- 
action because of the large number of harmonic components of 
the tidal potential and the reduced orbital separation near pe- 
riastron. By including the evolution of the MS star, especially 
of its rotation rate, many resonance crossings occur with en- 
hanced tidal interaction. We analyse the phenomenon of reso- 
nance locking whereby a particular tidal harmonic is kept res- 
onant with a stellar oscillation mode by the combined action of 
stellar evolution and other tidal harmonics. Resonance locking 
of prograde g-modes appears an effective mechanism for or- 
bital circularization of eccentric orbits. We consider the orbital 
evolution of the binary pulsar PSR J0045-7319 and conclude 
that resonance locking could explain the observed short orbital 
decay time of this system if the B-star spins in the direction 
counter to the orbital motion. 



of the components, the so-called equilibrium tide. However, in 
general the tidal deformation is time dependent and requires 



dynamical effects to be taken into account. Cowling (1941) 
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1. Introduction 

Two (or more) celestial bodies which orbit each other will 
experience periodic deformations due to the fluctuations of 
the mutually attracting gravitational fields. The tidally induced 
flows in these bodies are subject to dissipative processes, caus- 
ing a phase lag to develop between the perturbing potential and 
the displaced material. The resulting torque changes the orbital 
parameters of the system, and usually acts to drive the system 
towards circularization of the orbit and corotation of the com- 
ponents. If a binary system is in such a tidally relaxed state, it 
is relatively simple to calculate the tidally induced static shape 



noted that the tidal forces can excite gravity (g-) modes in the 
radiative layers of binary stars, the so-called dynamical tides. 
In case of a high mass early type star the stellar envelope is ra- 
diative and gravity modes can be excited for which significant 
damping occurs in a layer beneath the stellar surface, where 
the local therma l time scale is comparable to the period of os- 
cillation. Zahn ( 1977 ) developed an asymptotic theory valid 
for low forcing frequencies where the radiative damping is se- 
vere and resonances with g-modes are smeared out. However, 
when the system is far from equilibrium, tidal forcing contains 
relatively high frequencies, so that resonances with relatively 
weakly damped low radial order g-modes are likely to occur 
Zahn's theory may predict tidal timescales orders of magnitude 
too small when the asymptotic low frequency conditions are 
not met. 



In a series of papers Savonije & Papaloizou (1983; 1984) 
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presented non-adiabatic linearised calculations for perturbed 
spherical stellar models, enabling the study of resonantly ex- 
cited g-modes at different stages of evolution of the star To 
include effects of stellar rotation on the oscillations, an implicit 
2D hydrodynamics code has been developed (Savonije and Pa- 
paloizou, 1997) in which centrifugal distortion is discarded, but 
which fully accommodates the arising Coriolis forces in the ro- 
tating star. Due to the Coriolis force the response of a rotat- 
ing star to a spherical harmonic perturbing force is itself not 
a simple spherical harmonic function, but may be numerically 
solved on a 2D (r, i9) grid. Using the 2D-code, the effects of 
rotation on the stellar (g-)modes can be examined, while also it 
is possible to study the stellar quasi-toroidal r-modes (see Pa- 
paloizou and Pringle, 1978) in detail. In a previous paper (Witte 
and Savonije, 1999) (Paper I) we utilised the above mentioned 
2D code to calculate the g- and r-mode spectra of a 10 M© stel- 
lar model with a central hydrogen abundance of forty percent 
rotating uniformly at speeds up to forty percent of breakup. In 
order to study effects of the star's evolution, we extend the data 
of Paper I with some calculations for a Xc = 0.2 stellar model, 
the results of which are listed in the appendix of this paper The 
strength of the r-mode resonances, as given by the area of the 
peaks, becomes larger in a more rapidly rotating star. For rapid 
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stellar rotation with fig = 0.4 Vic, the strongest r-modes have 
peak areas comparable to g-modes with approximately 7-10 
radial nodes. Additionally the long wavelength weakly damped 
r-mode resonances occur at quite low forcing frequencies in a 
region where the stellar g-modes are severely damped as a re- 
sult of their short wavelength. If, therefore, a star is forced with 
strong harmonics that lie in this region, resonant excitation of 
r-modes can be expected to dominate the orbital evolution of 
the system. 



The recent discovery (see p<^aspi et al., 1994| ) of the binary 
radio pulsar PSR J0045-7319, which contains a B-star and a 
neutron star in a highly eccentric orbit, and for which radio 
pulsar timing has lead to a fairly accurate measurement of or- 
bital precession and orbital decay (see also Sect. triggered a 
number of publications in the field of dynamical tides in wide 
eccentric binary systems (e.g. Kumar and Goodman, 1996 ; Lai, 



1997; Kumar and Quataert, 1998 ). These papers considered 
some aspects of nonlinear coupling of modes, effects of stellar 
rotation and enhanced dissipation due to differential rotation 
by means of various approximations. We will now investigate, 
by detailed calculations, the effects of resonances with stellar 
oscillation modes (in a uniformly rotating MS star) on the tidal 
evolution of eccentric massive binary systems and apply our 
results to PSR J0045-7319. 

2. Basic equations 

We consider a binary system consisting of a Ms — 10 Mq main 
sequence star and a Mp = 1.4 compact (NS) companion 
in an eccentric orbit with eccentricity e and orbital period Poi b- 
Such a system could be easily produced during the supernova 
explosion following a period of mass transfer from the initially 
more massive primary to the rejuvenated MS object (van den 
Heuvel and Heise, 1972). 

The energy and angular momentum (magnitude) of the ec- 
centric orbit is given by 



Porb = -GMpMs/2a and H^.b - J^'f^ " 

where a is the semi-major axis and ui = 27r/Porb the mean 
angular velocity of the stars in their elliptic orbit. For simplic- 
ity we assume the stellar spin angular momentum vector Hg to 
be aligned with the orbital angular momentum vector. In this 
paper, we will only consider secular tidal changes in the mag- 
nitude of the orbital energy and angular momentum, and ignore 
changes in the configuration of the orbit as given by precession 
of the axes (which vanishes as the axes are assumed parallel) 
or advance of periastron (apsidal motion). 

2.1. The tidal potential 

Due to the motion of the companion in its eccentric orbit, the 
10 Mq MS star is exposed to a changing external gravitational 
field. We can facilitate the analysis by subdividing the forcing 
potential into its harmonic components and evaluating the con- 
tribution of each term to the tidal process separately. Labelling 



the companion's coordinates relative to the star with a prime, its 
perturbing potential is expanded as the real part of (e.g. Morse 
and Feshbach, 1952) 
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where P"'(cos i?) denotes the associated Legendre polynomial 
and tra — 1 for m = and 2 for m > 0. We assume m > 
and characterise retrograde relative orbital motion by a negative 
forcing frequency a — a — m^L^ in the stellar frame. Fourier 
expansion of the time dependent terms gives 



(7) 
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where M ~ tot is the mean anomaly. The Fourier coefficients 
(often called 'Hansen coefficients' in this context) are given by 



27r 



The relative elliptic orbit can be expressed in the parameters e, 
p and (p' or in e, a and the eccentric anomaly E: 



— = a(l — ecosP). 
1 + e cos ip' 

The eccentric anomaly E is related to M and (p' through (e.g. 



Brouwer and Clemence, 1961) 



E-esinE = M 
and 
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By scaling the Hansen coefficients as 



(cos I 



V a'+i J ^"'(/ + m)!'' 
the companion's tidal potential can be expressed as 

00 I 
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c ™ cos{nLut — nup), 



(2) 



i.e. each tidal spherical harmonic component {I, m) can be ex- 
pressed as a series of harmonics n, representing circular orbits 
with angular frequency nu). In this paper we only consider the 
dominant tidal components with / = 2 and to = or m = 2. 
The m — 1 contribution vanishes for the aligned case stud- 
ied here. Modes with to = have no p dependence, therefore 
no angular momentum can be transferred via these oscillations, 
only energy exchange can take place. For these axisymmetric 



M.G. Witte & G.J. Savonije: Evolution of eccentric orbits 



3 




Fig. 1. Hansen coefficients hn^^^ 
tricities e = 0.2, 0.4, 0.6 and 0.8. 



witli I 



2 for eccen- 



modes, the coefficients c'™ = c'™„. For m = 2 the coeffi- 
cients with n < are negligible. In Fig. [l| we plot the Hansen 
coefficients for m = 2 as a function of n for a few different 
values of the eccentricity. It is apparent that as the eccentric- 
ity increases, the number of contributing harmonics increases 
greatly, and also that the magnitude of the largest coefficients 
(which enters the expression for the torque quadratically) in- 
creases by a significant factor due to the diminishing orbital 
separation at periastron. 

In our calculations we include rt-values starting at n = 1 
and going up to rimax = 2{ni + I + l)/por to ensure that the 
companion's potential is accurately reproduced. For m = we 
take c'™ to be the sum of the positive and negative n contri- 
bution. Here, fpa = sj (ll^yi is the ratio between the orbital 
frequency in periastron ojpor and the mean orbital frequency tj. 

2.2. Orbital evolution by the tidal exchange of energy and 
angular momentum 

Situated in an eccentric orbit with its companion, the 10 Mq 
star is simultaneously forced at many harmonic frequencies 
(T„ = nu) — rnfls (in the stellar frame), the magnitude of 
each term being proportional to c'™. The harmonic compo- 
nents with n < no — mint(ils/cj) correspond to negative 
forcing frequencies which means that the corresponding fre- 
quency in the inertial frame is smaller than the stellar rota- 
tion frequency. Hence in the stellar frame these components 
excite oscillations which run backwards. Excitation of such 
retrograde modes gives rise to spin-down of the star, while 
excitation of prograde modes {n > uq) has the opposite ef- 
fect. This opens the possibility of simultaneous tidal interac- 
tion with counteracting resonances with prograde (spin-up) and 
retrograde (spin-down) oscillation modes. By assuming the re- 
sponse of the star to the tidal forcing can be approximated by a 
linear treatment the problem simplifies considerably since we 
may apply the superposition principle and decompose the stel- 
lar response into independent harmonic components. We fur- 
ther assume the stellar response to each harmonic has reached 



a steady state in which the tidal excitation is balanced by radia- 
tive and viscous damping in the oscillating star. These approxi- 
mations may break down during resonance passages although a 
steady state is usually a fairly good approximation, see discus- 
sion in Sect. ||. The steady state assumption allows us to calcu- 
late the stellar response as a strictly periodic phenomenon, i.e. 
we are not forced to follow the stellar oscillations on a dynam- 
ical timescale, which, at the current level and for the duration 
required to study secular evolution, would not be possible with 
present day computer facilities. The steady state tidal torque 
can be determined by applying the implicit 2D code mentioned 
above whereby the periodic term can be factored out. In Paper I 
we have shown that for each harmonic {I, m, n) in the forcing 
potential the work done by the tide (per time unit) on the star 
and the associated rate of change of spin angular momentum 
Hs can be expressed as: 

^i™ = a„7:!™ and i/i™ = mTj™ 



(3) 



whereby the torque integral is defined as 

7:!™ = -7rc':" / Y Im(p'(r,?9)) 
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(4) 



where Im stands for imaginary part and (cos ■&) is the asso- 
ciated Legendre polynomial of index m and degree I. In steady 
state E^™ equals the energy dissipation rate due to the n*'' har- 
monic of the tidal (Z, m) forcing. For a given forcing frequency 
(T„ — nu (in the inertial frame) the tidal perturbation of the 
stellar mass density p'{r, -d) occurring in the above integral fol- 
lows from the tidal response calculations in Paper I, see current 
appendix. Conservation of energy and angular momentum then 
implies that the rate of change of orbital energy and angular 
momentum follows by adding up the stellar rates of change in 
response to each harmonic term in the tidal potential and then 
reversing the sign: 



i,m n 



and 



l.m n 



By expressing the orbital eccentricity e in terms of the stel- 
lar masses and the orbital energy E'pib and angular momentum 
Horh (e.g. I^andau andLifshitz, 1959 ) 



e = Jl 



2Mt, 



-Borb-fforb 



with A/tot ~ Mp + Ms, we can express the rate of change of 
the orbital eccentricity as 



(1 - e2)Krb - cjVT^H, 



orb 



(6) 



We use this expression for the rate of change of the eccentricity 
to prevent numerical roundoff errors from producing negative 
eccentricities or generating eccentricity in a circular orbit. Fi- 
nally, the rate of change of the semi-major axis follows from 

Ida f GMpMA -' ■ 
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By defining the eccentric orbit for given (fixed) stellar 
masses through a and e the tidal evolution of the orbit can be 
followed by numerically integrating Eqs. (H)-(0) whereby at 
every timestep the tidal exchange rates Eorh and Horh are de- 
termined through Eqs. (|])-(^. 

To investigate the resonant exchange of energy and angu- 
lar momentum in detail, we developed a routine that enables 
us to calculate this tidal evolution by interpolating the tidal 
torque integral (Q), as a function of forcing frequency (T„, stel- 
lar spin rate fig and evolutionary state of the MS star (expressed 
in terms of the core hydrogen abundance Xc) in the data pre- 
sented in Paper I and Tables Bl and B2 (see appendix for more 



details). 

2.3. Resonant exchange of energy and angular momentum 

We have seen that in an eccentric binary system the tidal forc- 
ing of the star by its orbiting companion occurs in a non- 
harmonic time dependent manner Decomposition of the forc- 
ing potential into harmonic components introduces a range of 
frequencies at which different harmonic modes of oscillation 
in the star are excited simultaneously. These forcing frequen- 
cies change as orbital evolution progresses due to the fact that 
the orbit may shrink or widen, and due to changes in the stellar 
rotation frequency. Moreover, as the system evolves, the stel- 
lar oscillation spectrum itself changes due to restructuring of 
the star in response to stellar evolution and due to spin-up or 
spin-down of the star caused by tidal effects and the changing 
stellar moment of inertia. These stellar changes result in shift- 
ing of the eigenfrequencies of the stellar modes on the nuclear 
timescale of the star (which becomes relatively short, with the 
possibility of rapid tidal evolution, near the end of core hydro- 



gen burning (Savonije and Papaloizou, 1984)). Hence, as the 



orbit evolves, the forcing frequencies shift through the stellar 
oscillation spectrum, picking up resonances in the star as eigen 
frequencies are crossed. During a resonance passage, the rate 
at which energy and angular momentum exchange takes place 
between the star and the orbit (i.e. the companion) can be many 
orders of magnitude larger than in the non-resonant case. The 
total amount of exchanged angular momentum is obtained by 
integrating the torque over time; since the speed at which the 
forcing frequencies shift through the stellar oscillation spec- 
trum will generally be proportional to the torque, a relatively 
large peak area in the torque-frequency graph does not guaran- 
tee a large relative contribution to the tidal evolution process. 
The angular momentum exchange by a resonance passage of 
the n*'' harmonic component crossing resonance peak k can be 
expressed as: 



(AH), 



d(g-„-CTo,fc) 



d(cr„ - cro,fc) (8) 



dt 



with 7^'™ the torque integral (Q). The rate at which a cer- 
tain harmonic component n moves through the resonance with 
mode k (at frequency (To_fc in the stellar frame) is a function of 
the rate of change of the mean orbital and stellar angular speed 



and the rate of stellar evolution (characterised by the rate of 
hydrogen depletion in the stellar core Xc): 



dt 



AljJ 



= n— m— — 

dt dt 



daQ^k dXc _ daQ^k <i^s 
dXc dt dn^ dt 



After separating the dynamical and the stellar evolution com- 
ponent we obtain 



d(cr„ 



dt 



CTo.fe) duj 1 
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dt 

dapM \ d/, 



AH, 
dt 



dXr 



dXc 



dXc 
dt 



(9) 



with Is the stellar moment of inertia. We note that, while the 
second term on the right of Eq. (^ is proportional to iJorb, the 
first term is, for given harmonic n, fully determined by £^orb- 
Of course the last term is determined by stellar evolution. 

2.4. Resonance locking 

When a harmonic term n of the tidal potential comes into res- 
onance with a stellar oscillation mode the resulting exchange 
of energy and angular momentum between the MS star and 
its companion usually leads to a rapid shift through the reso- 
nance, so that the resonance condition is quickly lost and the 
tidal exchange is rather limited. However, under special cir- 
cumstances, the resonance condition may be sustained for a 
relatively long period of time. A necessary condition for such 
locking on to a resonance is that the right hand side of Eq. (^ 
attains a small value during resonance passage while, in order 
to drive a significant orbital evolution, the acting torque should 
be large. Substituting (^ and introducing the orbital moment 
of inertia /orb 



Af^+M '-El' ^'^^ expressed as: 



d(g„ - ao,k) 
dt 



eXc 



Im 



(10) 



where I = m = 2 (the m = contribution is usually negligi- 
ble). The dynamical factor is defined as 



Cn 



3ry 

-^orb 



m 



while the stellar evolution factor is 

/ , daoh\ dis daa,k 



dXr 



dXc 



Typically ~ 0.5 < m for to = 2 locking. For the 

orbits we will consider /orb is typically two or three orders 
of magnitude larger than /g. For sufficiently large values of 
n the diagonal factor C — Cnn, which governs the frequency 
shift relative to the resonance peak due to the action of the n'^ 
harmonic itself, attains a small value. A harmonic term in the 
forcing potential which has a small self-shifting coefficient ( 
must be driven deep into a resonance before its influence on 
the frequency shifting becomes prominent. If its self-shifting 
tendency draws itself through the resonance, the small value 
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of ( results in a relatively slow resonance crossing with corre- 
spondingly enhanced exchange of energy and angular momen- 
tum. However, if the self-shift works in the direction opposite 
to the resonance crossing, the frequency shift (due to the com- 
bined action of the other harmonics and the effect of stellar 
evolution) towards the resonance peak may be halted, resulting 
in a locked condition. This locked condition, whereby the self- 
shift balances the effect of all other harmonics, may last for a 
long period, therefore the total effect on the orbit can be sig- 
nificant. The most efficient locking will occur when C — 0, i.e. 

for the tidal harmonic n ~ n^; ~ mint^J^g^, which depends 
on the system's orbital period. The locking condition depends 
on the distribution of the Hansen coefficients, i.e. on the orbital 
eccentricity and on the location of the stellar oscillation spec- 
trum relative to the locked mode (on where the strongest modes 
are situated in frequency space). 

Consider the common case where the n"^ harmonic ct„ ap- 
proaches a resonance ao,k from below, i.e. (t„ — ao^k < 
whereby C = Cnn is slightly negative. Assuming a MS star with 
prograde spin and n^; > uq — m int(ils/w), where uq denotes 
the division between retrograde and prograde harmonics, we 
can define 



1^ C 

no<j<n 



> 



< 



with either to = or m = 2. Unless specifically stated, the 
initial evolution state of the MS star corresponds to a central 
hydrogen abundance = 0.4. 

Before we go on to the discussion of our numerical results, 
let us first express the rate of change of the orbital eccentricity 
(^ due to the action of a single resonance, by substituting the 
corresponding expressions (|]) for and ij^™ as: 



de^ 
dt 



OJ 



Ear 



orb 



Im 



(12) 



where the factor in square brackets is always negative. It fol- 
lows that prograde harmonics {n > uq) would require n < 
to/VI — to yield a resonant increase of the orbital eccen- 
tricity. Unless the eccentricity is unrealistically large these two 
conditions are mutually exclusive, so that resonances with pro- 
grade modes tend to decrease the orbital eccentricity. On the 
other hand, retrograde modes which have n < uq and T^™ < 
require n > m/ \/\ — to enlarge the orbital eccentricity. De- 
pending on the stellar rotation rate these two conditions can 
possibly be fulfilled, so that resonances with retrograde to = 2 
modes may enlarge the orbital eccentricity. Given the large 
number of harmonics that appear in orbits with high eccentrici- 
ties (of the order 10^ for e ^ 0.8), the orbital evolution of such 
binary systems can be quite complicated. However, before we 
study the high eccentricity case, we will first consider the sim- 
ple case of an almost circular orbit to gain some insight in the 
effect of resonances. 



EC 



> 0. 



Note that all harmonics with j < n have C„j < 0, while all 
higher harmonics have C,nj > 0, while the torque integrals 
q-im ^ Q ^Qj. ^ ^ ^jjj positive for n > uq. Since C < 

(but almost zero) the n'^ harmonic tends to move away from 
the resonance at a-Q^k and must be driven into resonance by the 
combined action of the other harmonics (which will interact 
with other resonances) and the weak (positive) effect of stellar 
evolution: 

-C^J™ <Si + S2 + S3 + eX,. 

Locking requires that, when driven deep into resonance, the n*^ 
harmonic begins to dominate all other effects, so that its action 
can prevent resonance passage: 



Cmax(Tj™) >Si + S2 + S3 + eX,. 



(11) 



The criterion is easily modified for a retrograde stellar spin or 
for n > Uq. We will meet several cases of resonance locking in 
the numerical results discussed below. 



3. Results for moderately wide binary systems 

From now on we express all frequencies in units of ilc and 
note that we consider only the strongest tidal components of 
the companion's potential, i.e. the component for which I — 2 



3.1. Case of a narrow, almost circular orbit 

First, we illustrate how the inclusion of resonant effects influ- 
ences the tidal evolution in case of an almost circular orbit, 
so that we have to deal with only few harmonics. This we do 
by evolving a fairly narrow orbit (Poib = 3 days) with ec- 
centricity e ~ 0.02 and a slowly rotating stellar component 
(Ps = 25 days) until corotation and circularization is estab- 
lished, see Fig. |[ In Fig. ^ we depict the initial location (fre- 
quency) of the stellar / = to = 2 modes in terms of the or- 
der n of the forcing harmonic and indicate which harmonics 
have the largest Hansen coefficients for this orbit. The to = 
modes play generally no prominent role in the evolution. For 
these low eccentricities (e^ ^ 1), we can use a low order ex- 
pansion of the potential of the orbiting companion in powers of 
the eccentricity to analyse our numerical results. To first order 



in e the Hansen coefficients can be expressed as h-^ 



3,0 



h{ = - 

potential 



ie, h. 
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1 and h 



3,2 



3e, 



1 



GMp /r\2 



4a \a 



jc, yielding for the tidal 



6eP2(cosz9) cos ujt 



(cos cos(wt — 2ip) 
+P|(cosi9) cos(2wt - 2ip) 
+ ^eP| (cos i9)cos(3wt - 2tp) ] + O(e^). 



(13) 
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Fig. 2. Orbital evolution of a weakly eccentric (e = 0.02) binary system which initially, at t = 0, harbours a slowly rotating 
MS star (Pg — 25 days) in an orbit with a period of three days. During tidal evolution, resonant interaction through the n = 3 
harmonic decreases the eccentricity to very low values (e = 0.00012 at t = 1.5 Myr), while non-resonant interaction via the 
n — 2 component forces the star into corotation. a Orbital (solid, left ordinate) and stellar (dotted, right ordinate) rotation period, 
b ratio of periastron distance to stellar radius, c eccentricity (solid, left ordinate) and ratio of stellar rotation frequency to orbital 
frequency at periastron (dotted, right ordinate), d timescale for the change of the eccentricity, e and f schematic representation 
of the frequency distribution of forcing harmonics and stellar resonance frequencies att — (panel e) and t = t,„ax (panel f). 
Here, frequency increases from left to right, the contributing harmonic potentials have been labelled n = 1 to 10 (= rimax) and 
the gray level of each cell centred on n reflects the magnitude of the corresponding Hansen coefficient h^'^. Due to the very low 
eccentricities only the cell with n — 2 (and for t — 0, very vaguely, n = 3) is printed gray. The frequencies of the stellar g-modes 
are given by thick lines, while the location of the r-modes is in between the dotted lines. In panel e, the r-modes are outside the 
range of the forcing harmonics and do not show up. The dashed line gives the frequency for which ( = (nn — (see Eq. (|lo|)). 
The white line (panel f) indicates the position where ct = 0. 
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Note that the sign of the Hansen coefficients is irrelevant since 
the tidal torque integral is proportional to the square of the co- 
efficients. The (here prograde) n — 2 harmonic is by far the 
strongest for low eccentricities. During orbital evolution, its ex- 
citation frequency (J2 only crosses the resonance frequency of 
the weak prograde gjQ-mode (i.e. the g-mode with / = 2 which 
has 20 radial nodes). This happens at t 7.2 x 10"' yr, only 
14 kyr after the (intrinsically weaker) n — 3 harmonic crosses 
the much stronger g^g resonance peak. Around this time, the 
orbital period shows a small dip, mainly due to the n = 3 reso- 
nance crossing. Because the g^o mode is the lowest frequency 
g-mode considered here, no further resonance crossings occur 
with the n = 2 harmonic during the remaining orbital evolu- 
tion. As the eccentricity drops during the orbital evolution, the 
n = 3 Hansen coefficient is further decreased and the follow- 
ing n = 3 resonance passages induce little or no features in 
the curve of the orbital period. Due to the off -resonant n — 2 
forcing frequency the star steadily spins up towards corotation. 

For the n'^ harmonic the time evolution of the eccentric- 



ity ( 12 ) can be expressed as 



de m — n 
= 2 

at Horh 



(14) 



Note that, initially, the star spins very sub-synchronously, so 
that the torque acting on the star mT^\™ > for all n. The 
n = 2 harmonic of the m = 2 potential does not give rise 
to significant changes in the eccentricity of the orbit because 
its circularizing effect vanishes to order O(e^). It appears that 
the potential term with n = 3 is responsible for the strongest 
circularization effect, while the n = 1 component could po- 
tentially increase the orbital eccentricity of the system. How- 
ever, whereas the frequency ai is so small that it does not 
cross any prograde resonance during the evolution of the or- 
bit (although it does cross the r-modes around t = 0.8 Myr), 
the n = 3 forcing frequency starts between the g-modes with 
fc = 12 and fc = 13, and passes through all of the modes 
§13 '^o §21^ the evolution progresses. Due to the small ratio 
{h^'"^ /h\' )^ = (|e)^ the total torque on the star is dominated 
by the off -resonant n ~ 2 contribution even when the n — 'i 
forcing frequency is close to an eigen-frequency of the star 
Therefore the rate at which the n — ?, frequency shifts 
through the spectrum varies only little. The g-modes with thir- 
teen to twenty radial nodes roughly have Act k, IO^'^ctq- As 
only a few percent of the resonance peak area for such a peak 
is further than 10 peak widths away from the peak location, 
while the peaks are typically 10^ peak widths apart, the tidal 
exchange rate due to the n = 3 harmonic, and therefore the rate 
of change of the orbital eccentricity, is not equally distributed 
in time. Indeed, from Fig. ^ it is seen that only during a res- 
onance passage the eccentricity of the system is significantly 
reduced, while only small eccentricity changes occur when CT3 
is in between resonances. 



3.2. Example of resonance locking in a moderately eccentric 
(e 0.25) 5 day orbit 

Let us now consider a moderately eccentric binary system 
with initially e — 0.25 and orbital period Poib = 5 days, 
which has a periastron separation of Dpa/Rs ~ 3.8. The ini- 
tial rotation speed of the MS star is adopted to be 0.8 times 
the angular speed of the orbiting companion in periastron, or 
Ps ~ 3.6 days. For this moderately eccentric orbit the har- 
monics of the tidal forcing run from n = 1 to n,„ax = 17 in 

= 
for 



n 



2: h. 



0.85, while /i„ < i /i- 



our calculation. The maximum Hansen coefficient at t 
corresponds to 

n 

is /s = 1.6 X 10^^ 



> 4. The orbital moment of inertia of the five day orbit is 
b = 9.1 X 10^^ g cm^, while the stellar moment of inertia 



only, the factor ^ 



g cm' 

la 



b/57. Since ^ 

2 / 9 dao.13 



-0.15 



which appears 



in Eq. ( p^ becomes positive for n = 9, so that we expect lock- 
ing to occur (if it occurs at all) with the harmonic n = 7 or 8, 
which appears indeed to be the case. Because ct„ and there- 
fore 7^'" is negative for n < uq — 2 and positive for n > 2, 
the locking conditions ( pi] ) apply when we consider locking of 
n = 7. Because of the relatively narrow and moderately ec- 
centric orbit only a few low-order harmonics play a significant 
role in the orbital evolution of this system and the sums Si in 
the locking conditions consist now of single terms. In Fig. |^ 
we plot the orbital evolution of this system. Three phases of 
rapid tidal dissipation can be discerned during which the or- 
bital period and eccentricity decrease substantially. During the 
first phase of rapid evolution, around i = 0.2 Myr, the initially 
sub-synchronously spinning star (fig = 0.8a;pcr) is spun up to 
almost pseudo-corotation, while the orbital eccentricity decays 
from e — 0.25 to w 0.19, in a few 10^ years. This is a clear 
sign of resonance locking. 

However, at t = no mode is resonantly excited, and the 
forcing frequency CT7, which is at the low frequency side of the 
prograde g^g mode, drifts towards the fc = 13 resonance peak 
mainly due to the combined effect of the star's evolution and 
driving by n = 2. This n — 2 harmonic excites a non-resonant, 
heavily damped low frequency (short wavelength) response at 
a frequency in between retrograde gjo and the r-modes. The 
n = 2 forcing thus induces a shift Act7 which is insensitive to 
the exact value of CT2 and, for a given magnitude of its Hansen 
coefficient, we can assume its effect to be approximately con- 
stant. Although the retrograde n — 2 forcing term gives rise 
to the strongest /re^Mency shift, its widening effect on the orbit 
is small compared to the orbital decay caused by the prograde 
forcing with n > 3, especially by the n — 7 harmonic which 
approaches the resonance with g^g. This relatively strong mode 
is forced with high n and thus efficient in bringing down the 
orbital eccentricity and orbital period. In Fig. || we zoom in on 
the gf 3 resonance passage of the n — 7 harmonic and plot the 
dominant contributions to the frequency shift of this harmonic. 

At about t = 0.1 Myr, CT7 runs into the g^3 resonance 
peak, and due to the negative sign of (^7,7 the associated tidal 
torque tends to shift the n = 7 harmonic in the negative di- 
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Fig. 3. Orbital evolution of a system with initial period Porb = 5 days and eccentricity e = 0.25. Initially the star rotates sub- 
synchronously (fig/wper = 0.8), but resonance locking of the n = 7 harmonic of the tidal forcing quickly forces the star to spin 
near (pseudo) corotation. During the subsequent tidal evolution of the system, a few more cases of resonance locking (now with 
the n = 6 harmonic) occur, recognisable by the more or less square shape of the dips in panel d, which shows the timescale for 
the orbital decay, a Orbital and stellar rotation period, b ratio of periastron distance to stellar radius, c eccentricity and ratio of 
stellar rotation frequency to orbital frequency at periastron, d timescale for the change of the size of the orbit, e and f schematic 
representation of the frequency distribution of forcing harmonics and stellar resonance frequencies at f = (panel e) and t = tmax 
(panel f). 
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Fig. 4. Locking of the tidal harmonic n — 7 onto the resonance 
with g^g. We plot the dominant contributions to the (rate of) 
frequency shifting of the n ~ 7 harmonic relative to the res- 
onance peak. The solid line gives the net rate of shifting by 
summing over all harmonics and adding the almost constant 
positive contribution of ~ 10^^^ s~^ due to stellar evolution. 
The dotted line represents A = (0-7 — (7o,i3)/Af7, i.e. the dis- 
tance to the peak of the resonance in units of the FWHM of the 
resonance peak. 

rection, back out of the resonance. However, the thermonucle- 
arly driven expansion of the MS star counteracts this tendency, 
while the retrograde n ~ 2 harmonic pushes the n = 7 har- 
monic even stronger into resonance. Although the gl^ reso- 
nance gives an intrinsically much stronger response than the 
non-resonant n ~ 2 response, the high order n — 7 forcing 
is significantly weaker ((/i2/^7)^ — 3.4 x 10^) so that this ef- 
fect is moderated. The n — 7 harmonic thus gets trapped by the 
combined action of stellar evolution and the retrograde n — 2 
harmonic which drive the system into resonance against the re- 
sisting self-shift of the locked harmonic. Thereby the n = 7 
harmonic remains in close resonance with g^g during a few 
times 10^ years, see dotted line in Fig. ^ Thanks to the rel- 
atively small value of jC?,?!, this locking occurs for a relatively 
large torque. Associated with the continuous resonant tidal in- 
teraction we see a rapid decline of the orbital eccentricity. This, 
however, weakens the tidal forcing of the n = 7 harmonic end 
strengthens that of the n — 2 harmonic. Thereby the n = 7 
harmonic is driven deeper and deeper into resonance with gf 3 
until it can no longer withstand the growing positive frequency 
shift due the n — 2 component and passes the top of the res- 
onance. Then the locking is broken and the n — 7 harmonic 
suddenly moves through and away from the resonance. A simi- 
lar, but shorter, locking phase occurs att = 1 .067 Myr between 
n — 7 and g^2- The dominant n = 2 harmonic crosses the r- 
modes at t — 1.375 Myr which can be observed in Fig. |^ as 
a single deep spike. Around t — 1.9 Myr, the orbit now be- 
ing narrower and C = for n ~7-8, a new prolonged locking 
phase sets in with n = 6 and gf^. Fig. |^ shows that signif- 
icant decay (~ 20%) of the orbital eccentricity occurs when 
a harmonic locks on to a (not particularly strong) resonance. 



During these relatively short (of order ~ 10^ years) locking 
phases the sub-synchronous star spins up considerably. How- 
ever, the retrograde n — 2 forcing, in combination with the 
continuous stellar expansion of the MS star brings the rotation 
rate down to small values. This continues until, near the end of 
the calculation, the stellar spin rate fig drops below the mean 
orbital speed uj after which the m = n = 2 harmonic becomes 
prograde. This brings the stellar spin down almost to a halt. 

3.3. Tidal relaxation of a 10 day orbit with high (e = 0.7) 
eccentricity. 

Let us now illustrate the enhanced effectiveness of the tidal in- 
teraction in a highly eccentric orbit by showing that a system 
which has a 10 day orbit and e = 0.7 will evolve towards tidal 
relaxation on a timescale short compared to the MS lifetime 
of the stellar component (which is approximately 2 x 10^ yr). 
Fig. H shows the evolution of such an orbit; at i = we adopt a 
stellar rotation rate equal to 0.7 times the orbital angular speed 
at periastron. As is seen from panel strong forcing occurs in 
the r-mode range of the stellar oscillation spectrum. Many pro- 
grade g-modes, including fairly strong ones, are in the wing of 
the Hansen distribution, creating many opportunities to trans- 
fer energy and angular momentum from the orbit to the star. 
Because of the location of = 0, far from any resonance, no 
locking is expected to occur during the initial stages of the or- 
bital evolution. 

Due to strong interaction with the stellar prograde g-modes 
the orbital period is halved during the first one million years. 
During this period the harmonics n = 9 to 5 cross the r-modes, 
keeping the stellar rotation rate subsynchronous. 

Between t = 1 Myr and t ~ 3 Myr the harmonic for 
which ~ shifts further away from the r-modes, making 
the remaining r-mode crossings less effective, into the region 
of the stellar oscillation spectrum where strong g-modes are to 
be found. As a result, resonance locking with prograde g-modes 
occurs, spinning the star up to supersynchronous speed. During 
this period the dominant retrograde forcing of the n = 2 har- 
monic, combined with the spin-down effect due to the stellar 
expansion, tends to slow down the stellar rotation. However, as 
long as resonance lockings of harmonics with higher n-values 
occur, the star is kept at a supersynchronous rotation speed. 
When around i = 3.1 Myr the m — n — 2 forcing frequency 
becomes prograde, the higher harmonics are no longer driven 
through the stellar oscillation spectrum in a direction opposite 
to their own shifting tendency, and equilibrium between two or 
more harmonics that produces resonance locking is no longer 
possible. At the end of the calculated evolution, when the sys- 
tem is significantly closer to the tidal equilibrium situation, the 
stellar spin-down due to stellar expansion and the spin-up due 
to the dominant low-frequency n = 2 forcing reaches approx- 
imate equilibrium, whereby the shift of the forcing frequen- 
cies through the stellar oscillation spectrum progresses on the 
(long) nuclear timescale of the MS star. We expect full relax- 
ation when the MS star approaches the end of core hydrogen 



10 



M.G. Witte & G.J. Savonije: Evolution of eccentric orbits 




1 2 3 

time (Myr) 



time (Myr) 




1 2 

time (Myr) 

rgl6-2() r-iTiodes C t *-* PS^"^ 



time (Myr) 



Pgl 




40 



50 



60 



70 



79 



Pg4 



10 



11 



Fig. 5. Orbital relaxation in a fairly wide orbit with high eccentricity. At t = the orbital period is Porb = 10 days, the 
eccentricity is e = 0.7, and the star has central hydrogen abundance Xc = 0.40. The star is initially rotating at seventy percent 
of the orbital angular speed at periastron. Within 3 Myr the eccentricity is reduced to ~ 0.08, and the orbit is shrunk to ~ 4 days. 
At i ~ 3.1 Myr the n = 2 component in the forcing potential becomes prograde with respect to the stellar rotation. When 
its off-resonant spin-up effect comes into equilibrium with the spin-down effect of the stellar expansion, the orbital evolution 
progresses on a nuclear timescale. At this time, the stellar rotation speed relaxed to ~ 0.85 times quasi-synchronisation. a Orbital 
and stellar rotation period, b ratio of periastron distance to stellar radius, c eccentricity and ratio of stellar rotation frequency 
to orbital frequency at periastron, d timescale for the change of the size of the orbit, e and f schematic representation of the 
frequency distribution of forcing harmonics and stellar resonance frequencies at t = (panel e) and t = tmax (panel f). 
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burning and the stellar evolution speeds up considerably see 
(Savonije and Papaloizou, 1984). 



4. Results for wide binary systems: application to the 
binary pulsar PSR J0045-7319 

Let us now consider systems similar to the SMC binary radio 
pulsar PSR J0045-7319, which has an orbital period of 51.17 
days and an eccentricity of 0.8080. This binary system con- 
tains a B-star which, for an assumed 1.4 Mq neutron star, has 
a mass (8.8 ± 1.8)M0. From the observed rather strong spin- 
orbit coupling the MS star is thought to be significantly de- 
formed by rotation, i.e. it is thought to rotate rapidly with its 
spin axis inclined with respect to the orbital plane. The peri- 
astron distance for such a system is only ~ 4 stellar radii, and 
pulsar timing revealed that the orbit decays on a timescale of 
^ ~ -S xlO'^yrCseelBelletal., 1995HLai et al., 1995|;K aspi 
et al., 1996). Therefore, this system constitutes an ideal test lab- 
oratory for the stud y of stellar ti des and neutron star birth kicks. 
According to Lai (1996; 1997), the short timescale of orbital 
decay suggests a significant retrograde stellar rotation with re- 
spect to the orbit. Indeed, Brandt & Podsiadlowski (1995 ) show 
from Monte Carlo simulations that after a supernova explosion, 
the spins of most stars in massive systems have large inclina- 
tions with respect to their orbital axes, and a significant frac- 
tion of systems (^ 20 percent) contain stars with retrograde 



spins. Kumar & Quataert (1997; 1998) argued that significant 



differential rotation, with a nearly synchronously rotating sur- 
face layer around a rapidly spinning interior, may be required 
to explain the observed timescale. 

Although we cannot, at present, take into account the in- 
clined orbit of PSR J0045-7319 we can study some aspects 
of this complex system. We will calculate the evolution of an 
aligned system with approximately the same orbital period and 
eccentricity, while adopting various values for the initial spin 
rate of the MS star: slightly super-synchronous rotation, highly 
super-synchronous rotation and retrograde rotation of the MS 
star With 'synchronous' rotation we mean here of course a ro- 
tation rate equal to the orbital angular velocity of the compan- 
ion at periastron. The initial eccentricity of e ~ 0.8 implies that 
the periastron frequency is o^pcr — 15 uj. For such a highly ec- 
centric orbit the tidal forcing is decomposed into a large num- 
ber (of order 10^) of harmonics, so that resonance crossings 
become ubiquitous (e.g. see Fig. ||d). 

4.1. Case a: slightly super-synchronous MS star 

Let us start with a 51.17 day orbit, eccentricity e = 0.808 and 
a stellar rotation rate fls = 1.1 ojpci, i-C- the MS star is set to 
rotate 10% faster than the largest angular speed of the com- 
panion. This corresponds to an initial stellar rotation period 
Pg ~ 3 days. In Fig. ^ it can be seen that the retrograde g- 
modes (k < 20) fall outside the tidal window. The dominant 
harmonics lie between 9 < n < 43, all other harmonics have 
Hansen coefficients at least 50% weaker The r-modes all clus- 
ter around harmonic n = 29, near the Hansen-peak and are thus 



expected to be prominent in the following tidal evolution. The 
prograde g-modes (no — 35) from g2o to g4 can be excited, 
although they cover the range rt > 62 and thus lie in the weak 
wing of the Hansen curve, so that the resonances will be weak. 
The zero point of ( depends, for given to = 2, only weakly on 
k and corresponds to n = nc^ 40. Because there is no stellar 
mode which can be excited by this harmonic, we do not expect 
resonance locking to occur Even when the stellar spin rate is 
decreased to zero the harmonics near C = can excite only 
high radial order g-modes (k ~ 18) which seem too weak to 
induce efficient locking. 

Let us now turn to the numerical results depicted in Fig. ^ 
As expected, very soon the rt = 29 harmonic of the I = m = 2 
tidal forcing comes into resonance with the fundamental r- 
mode rQ. After running through the whole range of r-modes 
(we consider A: = to A; = 10 in our calculations, see Table B2 
in Paper 1) the same process repeats itself with the n = 28 
harmonic, etc. until finally, at t ~ 4.3 x 10^ yr, the n = 20 
harmonic runs through the r-mode range. Each time a harmonic 
runs through the r-mode range a new bump appears in the curve 
plotted in Figs. ^ and b because the r-modes transfer energy 
and angular momentum from the rapidly spinning star into the 
orbit in such a way that the semi-major axis and the eccen- 
tricity tend to increase. Eq. ( p^ shows that for the initial ec- 
centricity e = 0.81 the harmonics with 5 < ti < no = 35 
will increase the orbital eccentricity, and the r-modes fall in 
this range. The effect of the evolutionary spin down of the star 
and of these repeated r-mode passages is that the star rapidly 
spins down to a sub-synchronous speed. The repetitive cross- 
ings of r-modes can also be seen in Fig. ||d. In between r-mode 
crossings the system now and then crosses a prograde g-mode 
(n ~ 70-80 pass through gio to gi4) giving rise typically to 
orbital \P/P\ — 10^ years. Due to the stellar spin down the 
prograde g-modes shift to lower frequencies so that stronger 
harmonics can excite them. This effect, in combination with 
relatively strong interaction with to = g-modes, becomes 
dominant after about t ~ 2 Myr whereby the orbital separation 
and eccentricity begin to decay (slowly). After 4 million years, 
the stellar rotation rate is reduced to approximately eighty per- 
cent of the pseudo-corotation rate (Porb — 4.2 days) while 
the eccentricity, after an initial rise, is diminished by a small 
amount. We note that the effect of stellar expansion and despin- 
ning r-modes prevent the MS star from attaining approximate 
corotation in periastron, even while the periastron separation is 
quite small in this system. 

4.2. Case b: highly super-synchronous MS star 

In this case we adopt the orbital parameters of case a, but the 
star is initially set to rotate at a higher rate: fis/wpcr = 1-8. 
We still consider the case of stellar rotation in the same direc- 
tion as the orbital motion of the companion. Fig. ^ gives the 
resulting orbital evolution, with the addition of a dotted line in 
panel ^ which gives the tidal timescale in case only forcing 
of the n = 48 (to = 2) harmonic would be taken into ac- 
count. During some 1.1 Myr the tidal timescale, and therefore 
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Fig. 6. Orbital evolution of a very eccentric (e = 0.808) binary system with a 51.17 day orbital period. At t = the star rotates 
supersynchronously (fig/wper = 1-1), but by the continuous stellar expansion and resonant interaction with the (retrograde) 
r-modes the originally super-synchronous star is spun down to ~ 0.8 times periastron frequency (dotted line panel c). Each time 
that one of the tidal harmonics n = 29 to n = 20 runs through the r-modes the orbital period and eccentricity increase. This 
corresponds to the repetitive bumps in panels a and b and in deep repetitive spikes in panel d. After t ~ 2 Myr the prograde 
g-modes have shifted to lower frequencies (panel e) and become increasingly more prominent, so that the orbital separation 
and eccentricity slowly begin to decay, a Orbital and stellar rotation period, b ratio of periastron distance to stellar radius, 
c eccentricity and ratio of stellar rotation frequency to orbital frequency at periastron, d timescale for the change of the size of 
the orbit, e and f schematic representation of the frequency distribution of forcing harmonics and stellar resonance frequencies at 
t = (panel e) and t = tmax (panel f). 
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Fig. 7. Orbital evolution of a very eccentric binary system. At t = the system's parameters Porb = 51.17 days and e = 0.808 
are the same as for the orbit in Fig. ||, but the stellar rotation speed in this case is higher; 1.8 times the companion's angular speed 
at periastron. Very soon the n = 48 harmonic in the potential runs into the fundamental r-mode and is locked near the resonance 
frequency for a period of more than one million years. During this period the star is spun down from 1.8 times periastron speed 
to less than 1.3 times periastron speed, while the orbital period is increased to approximately 72 days. During the locking the 
orbital eccentricity increases to almost 0.85. a Orbital and stellar rotation period, b ratio of periastron distance to stellar radius, 
c eccentricity and ratio of stellar rotation frequency to orbital frequency at periastron, d timescale for the change of the size of 
the orbit, e and f schematic representation of the frequency distribution of forcing harmonics and stellar resonance frequencies at 
t — (panel e) and t = imax (panel f). 
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Fig. 8. Locking of the tidal harmonic n = 48 onto the res- 
onance with rQ. We plot the rate of frequency shifting of the 
n = 48 harmonic relative to the resonance peak. The contri- 
bution due to stellar evolution is plotted dot-dashed, the contri- 
bution of the n = 48 excitation itself is dashed. The solid line 
gives the net rate of shifting by summing over all harmonics 
and adding the contribution due to stellar evolution; the dotted 
line represents A = (0-48 — aofi)/ Aa, i.e. the distance to the 
peak of the resonance in units of the FWHM of the resonance 
peak. 

the acting torque, is dominated by the single contribution of 
this n = 48 harmonic, indicating resonance locking. 

Fig. ^ displays the rate of frequency shifting of the n = 48 
forcing frequency relative to the r;^ resonance frequency. On 
average, the net rate of shifting is kept close to zero due to can- 
cellation of the stellar evolution (positive) and n = 48 (nega- 
tive) contributions which are in equilibrium for some 1.1 Myr 
The resonance passages of intrinsically strong (low radial or- 
der) prograde g-modes which occur (see Figs. ^ and Qf) in- 
duce only minor changes to the orbit and to the locking con- 
dition, due to their small corresponding Hansen coefficients, 
while passages of relatively strong harmonics through the ret- 
rograde g-mode resonances are unimportant because the modes 
which are crossed have many radial nodes and are thus in- 
trinsically weak through strong damping. Although some of 
these resonance crossings do show up as spikes in the solid 
curve in Fig. ^ their duration is too short to significantly in- 
fluence the resonance locking; many more of these spikes are 
not plotted due to graphical limitations. Only the passage of 
n = 49 through rg around t — 0.7 Myr is able to briefly push 
the locked harmonic further from resonance. During evolution, 
C48,48 is positive, but its value decreases as the orbit becomes 
wider The ability of the locked mode to resist being pushed 
through resonance therefore weakens, and a larger torque is re- 
quired to maintain the equilibrium, causing the harmonic forc- 
ing to drift closer towards the resonance peak. Eventually, the 
resonance is crossed and the locking is ended, which happens 
around t — 1.2 Myr. After this happens the CT48 forcing fre- 
quency drifts through the higher radial order r-modes, which 
have lower peak values and are unable to restore the locking 



condition. The n = 47 harmonic runs into the stellar rj^ mode 
before the 71 = 48 frequency reaches the last r-modes, and from 
this moment on the timescale (panel ^) is no longer dominated 
by the action of the previously locked harmonic, but instead 
shows a sequence of dips below the dotted line which mark 
the r-mode resonance crossings of n = 47. Somewhat later 
(t ~ 1.75 Myr), the next harmonic with n = 46 reaches the 
r-modes, causing a similar sequence of dips in the timescale 
curve. 

4.3. Case c: a retrograde spinning MS star 

We will now consider the case where the star is spinning in the 
sense opposite to the orbital revolution of its compact compan- 
ion, at thirty percent of its surface break-up speed. Furthermore 
we will adopt an initial orbital period Poib = 55 days, an ec- 
centricity e — 0.81 and a core hydrogen abundance Xc = 0.21, 
thus a fairly evolved MS star. For this initial configuration the 
tidal window is completely filled up with retrograde (in stellar 
frame) g-modes, see Fig. The peak of the Hansen distribu- 
tion lies around harmonic n = 21, while the ^-factor vanishes 
for n ~ 32 which is near the resonance with g^iQ. Thus un- 
der favourable circumstances this resonance may get locked, 
as will indeed be the case. Following the system until the core 
hydrogen content has been diminished to twenty percent, the 
orbital/spin evolution proceeds as depicted in Fig. ^ 

At about t = 0.162 Myr the n = 29 harmonic approaches 
the resonance with g^^^^ee Fig. The factor C29,29, or ^29 
for short as defined in ([lO|), being small from the very begin- 
ning, approaches zero and changes sign at t ~ 0.167 Myr. 
From this moment on, the self-shifting direction of the n = 29 
excitation is counter to the direction it has been shifting in, and 
it stops itself from passing the g^iQ resonance frequency. Also, 
since (29 is so close to zero, the torque associated with the 
resonant interaction hardly leads to a shift of the forcing fre- 
quency CT29 and the locking can be very efficient. It depends on 
the action of the remaining harmonics whether the locking will 
last. The lower order harmonics with n < 29 tend to push the 
n = 29 forcing in the positive direction, i.e. deeper into reso- 
nance while the harmonics with n > 29 tend to shift it back 
away from the resonance. Although the latter harmonics can 
excite strong low radial order g-modes, the peak of the Hansen 
distribution occurs for n ~ 20, so that the low order harmonics 
dominate and keep the n = 29 forcing close to resonance, see 
Fig. |l^. The dotted curve in Fig. ^ and Fig. [l0| b represents the 
sole contribution of the n = 29 harmonic to the decay time of 
the orbital period. It can be seen that from t ~ 0.164 Myr to 
t ~ 0.171 Myr most of the orbital decay is caused by the reso- 
nant interaction between n = 29 and g^iQ. During some 1 3 kyr 
the orbital decay time is between 10^ and 5 x 10^ years, which 
is the observationally estimated decay time in PSR J0045-73 19. 
Later in the evolution the n = 29 harmonic causes more phases 
with a similar short decay time while exciting g^iQ, but the 
ctq.-io resonance frequency is never crossed. 

In our calculation with retrograde but aligned spin lock- 
ing occurs for conditions slightly off the observed Poib and 
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Fig. 9. Orbital evolution of a system mimicking PSR J0045-7319. Alt = the orbital period is Porb = 55 days, the eccentricity 
is e = 0.81, and the star has central hydrogen abundance = 0.21 so that the periastron distance is initially approximately 
3.8 stellar radii. The star initially has a retrograde rotation frequency of Og — 0.30c- The large number of contributing tidal 
harmonics in the early stages of the orbital evolution frequently cause the timescale for orbital evolution (panel d) to reach 
values below 10^ yr. A total of approximately 1000 resonance passages occur during the evolution shown, which is terminated 
as Xc reaches 0.2. a Orbital and stellar rotation period, b ratio of periastron distance to stellar radius, c eccentricity and ratio of 
stellar rotation frequency to orbital frequency at periastron, d timescale for the change of the size of the orbit, e and f schematic 
representation of the frequency distribution of forcing harmonics and stellar resonance frequencies at i = (panel e) and t = imax 
(panel f). 



16 M.G. Witte & G.J. Savonije: Evolution of eccentric orbits 




Fig. 10. Resonance locking of n = 29 forcing harmonic on the retrograde g^j^Q-mode. a Left ordinate: net rate of relative 
((T29 — CTo.io) frequency shifting due to the combined effect of all tidal harmonics and stellar evolution (solid line) and the 
shifting rate due to excitation of the n — 29 term alone (dashed line). The nature of the resonance crossings is indicated along the 
solid line as ng^k, where n gives the forcing harmonic and k the number of radial nodes of the excited oscillation. Right ordinate: 
(29 (dotted line) and the distance between the n = 29 forcing frequency and the fc = 10 retrograde g-mode resonance frequency, 
b Detail of panel characteristic timescale of orbital evolution (solid line) and characteristic timescale due to n = 29 forcing 
alone (dotted line). 



M.G. Witte & G.J. Savonije: Evolution of eccentric orbits 



17 



e. But note that the system parameters for which locking oc- 
curs depend on the details of the stellar input model and on the 
procedure for interpolating between the two evolution states 
Xc = 0.4 and Xc = 0.2. We conclude that resonance locking 
might be responsible for the observed fast orbital decay in this 
system. 

5. Discussion 

We have studied the tidal evolution of eccentric early type bi- 
nary systems by decomposing the time dependent tidal poten- 
tial in a series of harmonics and by calculating the combined 
linear response of the uniformly rotating MS star with help of 
2D-implicit code. We assumed the tidal oscillations to be in a 
steady state, whereby the tidal excitation is balanced by radia- 
tive damping. This allowed us to follow the orbital and spin 
evolution of the binary system in detail. However, during res- 
onance crossings the assumptions of steady state and linearity 
may break down. Resonance peaks may sometimes be crossed 
on timescales short compared to their corresponding damping 
timescale rj. As build-up of the full steady state response will 
take a number of damping times one may question the validity 
of the steady state approximation under these circumstances. 
In reality the resonance crossing will be slower, as the tidal re- 
sponse is not yet fully developed. However, the tidal evolution 
is not expected to be much different from our steady state re- 
sults because the real evolution timescale is stiU likely to be sig- 
nificantly shorter than the timescales associated with the other 
(non-resonant) harmonics. 

Non-linear effects can become significant during a reso- 
nance crossing with a weakly damped mode and may give rise 
either to enhanced or reduced tidal evolution, depending on cir- 
cumstances. But note that non-Unear effects may often be of 
less importance because strong tidal interaction generally im- 
plies rapid orbital and spin evolution away from the resonance, 
so that the oscillation amplitude is already limited by linear 
effects. Non-linear effects may influence resonance locking, 
whereby a harmonic remains near a resonance for a prolonged 
period of time. Non-linear damping might lower the resonance 
peak height and might thus diminish to some extent the dura- 
tion of locking. 

On the other hand, we have ignored the effect of resonances 
with I > 2. In a rotating star the dominant I = 2 tidal compo- 
nent can excite (through the Coriolis force) also I =4, 6, 8 etc. 
modes, see Paper I. These higher spherical degree resonances 
will add additional resonance crossings during tidal evolution 
although for moderate stellar spin rates the peak area of these 
resonances is much less than that of the I = 2 peaks. How- 
ever, since peak area is not important for resonance locking, 
resonance locking may be even more important than our calcu- 
lations show due to the denser oscillation spectrum. 

6. Conclusion 

Our calculations indicate that the many forcing harmonics 
present in orbits with significant eccentricity give rise to ef- 



fective tidal evolution. By taking into account not only the or- 
bital evolution, but also the evolution of the MS star - espe- 
cially of its rotation rate - the system progresses through many 
resonances with stellar oscillation modes, a process that drives 
itself. Frequently conditions are favourable for resonance lock- 
ing, whereby a particular tidal forcing harmonic is kept near- 
resonant with an oscillation mode for a prolonged period. Such 
phases lead to rapid orbital and spin evolution. It appears that 
fairly wide massive binary systems with large eccentricity can 
be almost circularized in a few million years. The observed 
short orbital decay time of the binary pulsar PSR J0045-7319 
could be due to resonance locking if the stellar spin is retro- 
grade. 
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I = 2 



a = 0.0 f^. 



Xr = 0.2 




Fig.Al. Torque integral T'™ versus forcing frequency a for 
forcing with I = 2 and m = on a 10 Mq non-rotating star 
Prograde g^ -modes with k — 4 and 11 are labelled. Crosses 
denote calculated points, while the solid curve represents a fit. 



Appendix A: Determination of the tidal torque 

In Paper I we calculated the non-adiabatic tidal response of 
a uniformly rotating, somewhat evolved (Xc — 0.4), 10 Mq 
main sequence star to the dominant I = 2 components of its 
compact 1 Mq binary companion's tidal potential for a wide 
range of prograde and retrograde forcing frequencies and de- 
termined the corresponding tidal torques on the MS star The 
relative orbit of the compact star was assumed circular and with 
a fixed separation of 4 stellar radii (of the MS star). In order 
to include effects of stellar evolution on the tidal evolution of 
the orbit, we have now calculated additional torque spectra for 
a stellar model with core hydrogen abundance (mass fraction) 



Xr — 0.2, the results of which are listed in Tables Bl and B2 



Both stellar models, comprising 1200 (radial) zones, were 
calculated using the stellar evolution code developed by Eggle - 
ton (1972) and more recently revised by Pols et al. (1995), 



with the addition of recent OPAL opacity Tables (Iglesias and 
Rogers, 1996). 

The model with Xc = 0.4 has a stellar radius equal to Rs = 
3.825 X 10" cm, an effective temperature T^s = 2.314 x 10^ K 
and a stellar moment of inertia Is ~ 1.55 x 10^^ g cm^. The 
break-up angular speed equals fl^ — 1.54 x 10^"* s^^. The 
model with X^ = 0.2, on the other hand, has a radius = 
4.807X 10" cm, an effective temperature T^s = 2.157x 10'' K, 
moment of inertia /g = 2.059 x 10^^ g cm^ and break-up an- 
gular frequency ilc — 1.095 X 10^^ s^^. 



Fig. Al shows the calculated torque integral spectrum for 



m = forcing on the non-rotating star; compared to the 
Xc = 0.4 data, the degree of central condensation in the star 
has increased, the resonance frequencies of the stellar g-modes 
shift to much higher dimensionless values CTo/f^c- The break- 
up frequency ilc however, drops due to the increasing stellar 
radius, such that the actual frequency of the modes on average 
increases by only a few (~ 5) percent. It is also seen that the 



spacing between the individual modes is very different for the 
two stellar models. Different modes react differently upon the 
changes in the stellar structure, resulting in seemingly random 
shifting of the peaks superimposed on the global trend towards 
higher dimensionless frequencies. 

For Xc = 0.2 we synthesise g-mode spectra for fig — 0.2, 
0.3 and 0.4 flc by normalising the peak widths and areas from 
the TO = 0, rig = 0.0 ftc, Xc — 0.4 spectrum from Paper I to 
the rn = 0, Jig = 0.0 57c spectrum for Xc = 0.2 and assuming 
that the deducted ratios apply to the other stellar rotation rates 
as well. We also assume that the coefficients c, for the il,- 



of the resonance frequencies 



dependence ^ 

calculated for the Xc = 0.4 model in Paper I can be used for 
the Xc = 0.2 spectra. 

At low forcing frequencies where the resonant g-mode os- 
cillations attain high radial orders the values of T'™ are dom- 
inated by the turbulent dissipation in convective regions (see 
Paper I). For the orbital evolution calculations of this paper we 
approximate the torque integral in these frequency regions by a 
small constant level. For the non rotating Xc = 0.4 model this 
level is 1 X 10^ erg, for the Xc = 0.2 model it is 5 x lO^'"^ erg. In 
a rotating model the calculated torque values in the retrograde 
frequency domain around the r-modes become larger Here we 
multiply the level with a factor (l + 10ris/57c) to approximately 
fit this effect. The area ActT'™ of this low frequency interval of 
the stellar spectrum is not very large, and compares to the peak 
area of a g-mode with approximately 15 radial nodes; much 
less than the peak area of the low radial order g-modes (1 to 
10 radial nodes). Under many circumstances, this regime will 
therefore not play an important role in the evolution of the or- 
bit. Under some circumstances however, e.g. if many strong 
tidal harmonics have frequencies in the designated area, or if 
only few harmonics are present in the forcing potential and the 
strongest one has its frequency in this domain (as is the case 
during resonance locking in moderately eccentric orbits, see 



Sect. 3.2), the details of our orbital calculations are subject to 



uncertainties due to the numerical limitation of having to work 
with a finite mesh and due to the theoretical limitation of not 
having a sound theory for turbulent convection. Nevertheless, 
the calculations presented here are believed to give a qualita- 
tive picture of the process of tidal evolution and are useful as 
an order of magnitude estimate for the strength of the occurring 
effects. 



A.l. Interpolation procedure for the tidal torque 

For each forcing frequency (t„, the corresponding torque inte- 
gral T^J™ is calculated by interpolation in the Tables of Ap- 
pendix B: and those in the appendix of Paper I in terms 
of the three state parameters: forcing frequency a-„, stel- 
lar rotation rate ilg and evolutionary state of the star Xc 
(i.e. the hydrogen content of the convective core). Subse- 
quently the obtained torque values are normalised by a factor 
(A/p/M0)2(4i?,/a)(2'+2)(/ii'+i)'™^2 



as the torque integrals 
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have been calculated for a 1 Mq companion in a fixed circu- 
lar orbit with a = AR^. 

To perform interpolations in fig, linear relations were used 
for A(T and the peak area ttActT^j'™, while polynomial in- 
terpolations of the highest possible order were used for the 
resonance frequencies a^. For a^, we assumed the deriva- 
tive ^f- for m = modes to vanish for a non-rotating star 
(from symmetry), and for r-modes to become — | in order to 
obey the asymptotic value a — —j^^^ls with / = 3 and 
m = 2. The r-mode resonance frequencies were fitted by 

n~ 3 



C2X^ + c^x^ + C4X^ + c^x^ with X = rU^, 



whereby the polynomial coefficients are listed in Table B3 



Linear interpolations are used between the X,, = 0.4 
and Xc = 0.2 models for the peak widths and areas and the 
frequency ctq and polynomial coefficients Cj of each mode n in 
the spectra, and for the stellar radius and moment of inertia. 

As the calculated resonance peaks are accurately repre- 
sented by the harmonic oscillator fits, the peak width of the 
fits can be assumed to give a good measure of the damping 
timescale ta = of the excited modes. From the tabulated 
peak widths we find that for all calculated g-modes Td lies be- 
tween ~ 45 days and ~ 10^ yr; for the r-modes Td varies be- 
tween ~ 200 days and ^ 200 years unless the stellar rota- 
tion drops below 0.1 fie in which case both the strength of the 
modes and the damping rate is strongly reduced. 

While evolving the orbit, timesteps are chosen in such a 
manner that no frequency component f7„ will travel more than 
what corresponds to five percent of the height of a resonance 
peak per timestep. 
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Appendix B: Tables 

Table Bl. Fitting parameters for g^-mode resonances with an ? = 2, to = tidal potential in a 10 Mq star with = 0.2. All 
tabulated values for the torque integral have been obtained for a fixed orbital separation. 



= 0.0 fie 



k 


fTo/Oc 


(erg) 






k 


(TO /He 


To"" (erg) 






k 


tro/Oc 


-71,"" (erg) 






1 


3.4254 


1.23 X 10""^ 


1.14 X 10" 


- 5 


8 


1.0127 


1.94 X 10*^ 


2.06 X 10" 


-4 


15 


0.5382 


2.63 X 10*" 


8.91 X 10" 


-4 


2 


2.4130 


4.53 X 10*^ 


1.16 X 10" 


-5 


9 


0.9356 


7.18 X 10*^ 


1.00 X 10" 


-4 


16 


0.5142 


2.45 X 10*° 


7.41 X 10" 


-4 


3 


2.1362 


5.93 X 10** 


2.93 X 10" 


-6 


10 


0.8057 


1.28 X 10*^ 


1.81 X 10" 


-4 


17 


0.4768 


2.30 X 10*° 


5.76 X 10" 


-4 


4 


1.5940 


5.29 X 10** 


2.34 X 10" 


-5 


11 


0.7592 


4.86 X 10*^ 


2.97 X 10" 


-4 


18 


0.4528 


1.21 X lO^''' 


1.91 X 10" 


-3 


5 


1.3853 


3.32 X 10*^ 


4.26 X 10" 


-4 


12 


0.6759 


7.13 X 10*^ 
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-4 


19 
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3.13 X 10^^ 
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-4 


6 
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-2 


13 
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-4 


20 
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-3 


7 


1.1867 
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-4 


14 


0.5846 


2.35 X 10*^ 


2.43 X 10" 
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Table B2. Fitting parameters for r|-mode resonances with an / = 2, to = 2 tidal potential in a 10 M© star with X^. = 0.2. All 
tabulated values for the torque integral have been obtained for a fixed orbital separation. 
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X 
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-2 
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-4 


-1.1721 


X 
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-1 
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X 
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-2 
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-4 
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X 
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-1 
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-3 


6 
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X 
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-4 
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X 
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-1 
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-3 


7 


-8.4114 


X 
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-2 
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-4 
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X 
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-1 
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8 
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X 
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X 


10" 


-2 
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-3 


9 
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X 
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-2 
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-4 
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X 
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-2 
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1.44 X 10" 


-3 


10 


-7.5465 


X 


10" 


-2 


-1.01 X 10*^ 
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-4 


-8.6386 


X 


10" 


-2 


-1.17 X 10^^ 


1.69 X 10" 


-3 
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Table B3. Coefficients Cj for the fig-dependence ctq = Cj {^jf^ of the resonance frequency of r-modes in a 10 Mq star with 
Xc = 0.4 and = 0.2. (ci = -i) 



Xe = 0.4 



k 


C2 






C3 




C4 




C5 











9.76 X 


10" 


-3 


-3.29 X 10" 


-2 


3.61 X 10" 


-1 


-4.23 


X 


10" 


-1 


1 


1.39 X 


10" 


-2 


-3.55 X 10" 


-2 


4.32 X 10" 


-1 


-4.42 


X 


10" 


- 1 


2 


9.65 X 


10" 


-3 


7.34 X 10" 




3.53 X 10" 




-5.10 


X 
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3 


3.69 X 


10" 


-3 


2.73 X 10" 


-1 


-5.79 X 10" 


-2 


-1.51 


X 


10" 


- 1 


4 


2.25 X 


10" 


-4 


4.77 X 10" 


- 1 


-4.43 X 10" 
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X 


10" 
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5 
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10" 


-3 


6.25 X 10" 
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-2.19 


X 


10" 




6 


-1.38 X 


10" 


-3 


8.98 X 10" 




-8.40 X 10" 




-1.22 


X 


10" 




7 


-5.59 X 


10" 


-3 


1.23 




-1.59 




3.45 


X 


10" 


- 1 


8 


-1.27 X 


10" 




1.64 




-2.73 




1.33 








9 


-1.60 X 


10" 




2.00 




-3.68 




1.99 








10 


-3.45 X 


10" 




2.68 




-6.37 
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